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We introduce the possibility to use characteristic natural neutrino backgrounds, such as Geo- 
neutrinos (p e ) or solar neutrinos (v e ), with known spectral shape for the energy calibration of future 
neutrino detectors, e.g. Large Liquid Scintillator Detectors. This "CalEffect" could be used without 
the need to apply any modifications to the experiment in all situations where one has a suitable 
background with sufficient statistics. After deriving the effect analytically using \ 2 statistics, we 
show that it is only tiny for reactor neutrino experiments, but can be applicable in other situations. 
As an example, we present its impact on the identification of the wiggles in the power spectrum of 
supernova neutrinos caused by Earth matter effects. The Self-Calibration Effect could be used for 
cross checking other calibration methods and to resolve systematical effects in the primary neutrino 
interaction processes, in particular in the low energy cross sections. 



1. INTRODUCTION 

As neutrino physics is entering the stage of precision 
measurements, future detectors will have to face new 
challenges in terms of background reduction, control of 
systematical uncertainties, and calibration of the count- 
ing efficiencies and the energy reconstruction. The most 
widely used techniques for performing the energy cal- 
ibration are radioactive sources with known properties, 
artificially accelerated lepton beams, or indirect methods 
based on the interaction products of neutrinos (see e.g. 
Refs. [U, H, El ) • All of these methods have in common that 
they arc only sensitive to charged particles, i.e. possible 
uncertainties associated with the primary neutrino inter- 
action that produces these particles are not taken into 
account. Therefore it is desirable to calibrate the detec- 
tor directly with a neutrino beam. Such a measurement 
has been performed by GALLEX using a Cr-51 neutrino 
source, but it was very involved and required a special 
detector design [1] . The SAGE experiment has also been 
calibrated with neutrinos, from Cr-51 @ and Ar-37 @|. 

In this letter, we propose to use natural neutrino 
sources with known characteristic spectra, in particular 
Geo-neutrinos and solar neutrinos as calibration sources 
for future low-energy neutrino detectors. These types of 
neutrinos are always present, and, in a several-tens-of- 
kilotons detector, they induce at least a few thousands 
of events per year. The energy spectra for both, Geo- 
neutrinos (z7 e ) and solar neutrinos (f e ), are known very 
precisely since their production processes are very well 
understood. They exhibit characteristic steps at the cut- 
off energies of various processes, which are easy to locate 
in the final event spectrum even if the total rate is not 
known. 



a Email: Joachim.Kopp@mpi-hd.mpg.de 
b Email: Manfred.Lindner@mpi-hd.mpg.dc 
c Email: Alexander. Merle<Qmpi- hd.mpg.de 



Of course, the energy calibration with Geo-neutrinos 
or solar neutrinos will not eliminate the need for other 
calibration methods, in particular because it is insensitive 
to spatial or temporal variations of the detector perfor- 
mance and can be used only in the very low energy re- 
gion, but it is nevertheless interesting since it constitutes 
a self- calibration of the detector and provides an inde- 
pendent cross check of other methods without requiring 
any modifications of the detector design. 

In Sec. [H we present the small effect of the detector 
self-calibration on reactor neutrino experiments and ex- 
plain the effect analytically. Then, in Sec.[3l we will show 
how the Self-Calibration Effect can improve the accuracy 
of supernova neutrino measurements, and we finally con- 
clude in Sec. 0] 

2. BACKGROUND NEUTRINOS AS 
CALIBRATION SOURCES 

A future Large Liquid Scintillator Detector (LLSD) 
will have very good statistics even for background neu- 
trino sources such as distant nuclear reactors and Geo- 
neutrinos. This opens the interesting possibility to use 
these neutrinos as calibration sources for the energy re- 
construction, which leads to the seemingly paradoxial sit- 
uation that a measurement with backgrounds can yield 
more precise results than a measurement without back- 
grounds. 

Geo-neutrinos are particularly well suited for this self- 
calibration because their spectrum, which is given by 
the uranium and thorium decay chains, has character- 
istic steps at very well known energies that can easily be 
located in the data samples (see Fig. [!} . An analogous 
self-calibration for i> e 's may be possible using solar neu- 
trinos since the spectrum of Be-8 and pep is also very 
characteristic. 

Of course, background neutrinos are insensitive to 
short-term variations of the detector properties such as 
temperature fluctuations, so they will not eliminate the 
need for other calibration methods. But they are the 
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only tool to resolve systematical effects in the primary 
neutrino interaction process such as uncertainties in the 
low energy neutrino-nucleon cross sections. 

In this chapter, we will first decribe in which analysis 
we have found the self-calibration effect and then show 
that it can be explained analytically using the x a P~ 
proach. 

2.1. The tiny Self-Calibration Effect in Reactor 
Neutrino Experiments 

Already in 0] we have discussed the physics potential 
of reactor experiments with a LLSD such as the proposed 
45 kt LENA detector @, Q. Using the GLoBES soft- 
ware [IU , we have investigated the accuracy of such 
an experiment to measurements of the reactor angle #13 
as well as on the solar oscillation parameters, #12 and 
Amy, using mobile nuclear reactors amongst others. As 
backgrounds, we have taken into account the 20 closest 
reactors to the possible LENA site Pyhasalmi in Finland 
as well as Geo-neutrinos coming from uranium and from 
thorium (the Geo-neutrinos from potassium are not in- 
cluded in the analysis, since their energies are below the 
threshold for inverse beta decay, which is used as detec- 
tion reaction for the reactor-zVs in a LLSD, see Fig.[T]). 
For 0i3, these backgrounds turn out not to be very impor- 
tant: the optimum baselines for such a measurement are 
very short (~ 1 km), which leads, due to the large fidu- 
cial mass of the considered detector, to extremely high 
event rates even for a mobile reactor. Hence, for a mea- 
surement of the small reactor angle, the perturbation by 
the two different backgrounds is so small that it is not 
possible to exploit any information on the spectral shape 
of Geo-neutrinos. 

However, for a measurement of #12, the situation is dif- 
ferent: considering the SMALL scenario (details on this 
analysis can be found in Q), which corresponds to a mo- 
bile nuclear reactor with a thermal power of 0.5 GWth 
and 2 years of data taking, one can see a small effect 
induced by the background self-calibration of the experi- 
ment, which will be pointed out here. For the inclusion of 
Geo-neutrinos in our analysis, we have considered three 
different situations: 

• No Geo-neutrinos: In this case, Geo-neutrinos 
are completely absent, and only the background 
from distant nuclear reactors is taken into account, 
which yields 850 events per year. 

• Geo-neutrinos with a 10% uncertainty: Here, 
Geo-neutrinos coming from uranium and from tho- 
rium are taken into account and the uncertainty 
in both their fluxes is assumed to be 10%. The 
Geo-neutrino spectra in our simulations are taken 
from [H, [l| . The reactor background is the same 
as in the scenario without Geo-neutrinos. 

• Geo-neutrinos with a 100% uncertainty: This 
scenario is equivalent to the previous one, but now 
the uncertainties in the two Geo-neutrino contribu- 
tions are taken to be 100%. 



The result of our analysis for the SMALL scenario is 
again plotted here in Fig. [2] (cf. Fig. 3 in 0]). The 
interesting region is the marked rectangle. Taking a 
closer look it turns out to show a seemingly paradox- 
ial situation: from Fig. [3] one can see that, for certain 
baselines, the lower bound on the assumed true value of 
sin 2 2#i2 = 0.83 is slightly better with the background 
by Geo-neutrinos (and even in both cases, with 10% and 
100% flux uncertainty) than without. The reason is that 
this background has such a characteristic form due to 
the well-known Q-values of the nuclear production re- 
actions of the neutrinos that the positions of the steps 
in the spectrum can, in the minization of the \ 2 func- 
tion, be used to reduce a systematical bias in the energy 
calibration essentially to zero, which may not be possi- 
ble without such a characteristic background. Hence it 
can indeed be the case, that, for certain configurations, 
a suitable background can improve, rather than worsen, 
the sensitivity of certain kinds of experiments. 

The effect for reactor neutrinos is very tiny. This is 
because one needs to have large background rates to be 
able to resolve the characteristic spectrum properly. Also 
the energy scale of the background is crucial since it may 
not be possible to extrapolate a calibration made by MeV 
neutrinos accurate enough e.g. to the GeV scale. How- 
ever, the principle effect existis, as we will demonstrate 
in the next section (and what we have tested using dif- 
ferent hypothetical benchmark scenarios), and one just 
has to find a scenario where this effect can be applied. 

2.2. Analytical discussion of the Self-Calibration 
Effect 

Let us now discuss the self-calibration effect analyt- 
ically. We consider a simplified \ 2 expression using a 
Gaussian approximation for the distribution of the event 
rates: 

v 2 \ (Ti(a h b) - Nj) 2 qfl 

1 L _ 

Here, the Ti's are the theoretically predicted event rates 
calculated with certain test values for the neutrino os- 
cillation parameters and the systematical errors, while 
the A^j's are the rates obtained by using the assumed 
"true" values for these parameters. The a.j's are nui- 
sance parameters as e.g. overall normalization errors or 
bin-dependent shape errors. The explicit form of the T,'s 
is 

T^Q. + ajNtlb), (2) 
N i (b) = (l + b)-[(N m)i+1 -N 5(i) )- 

■(*(<) -L*(<)J)+%(i)jL (3) 
8(i) = b-{i + t {) + l -)+i. (4) 

Here, Ni(b) are the rates for wrong energy binning im- 
plied by a non-zero energy calibration b which are ob- 
tained from the correctly binned rates Ni according to 
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Eq. ([3]), which is essentially a linear interpolation be- 
tween the events in bin [S(i)\ + 1 and bin |£(*)J> is 
the energy threshold of the detector, expressed in terms 
of the bin width, and the Gaufi bracket ['J denotes the 
floor function. 

The normalization uncertainties are completely con- 
tained in the bin-to-bin normalization factors eti with 
errors <7j which contain all spectral distortion effects, 
such as imperfect knowledge of the reactor signal spec- 
trum, neutrino oscillations, etc.. The energy calibration 
is parametrized by b. Normally, one would also have to 
include a penalty term for the energy calibration of the 
form b 2 ja\ (as done for the spectral distortions). Such a 
term is omitted here because we assume no external in- 
formation on 6, which means that the energy calibration 
is initially completely arbitrary. We do, however, assume 
b to be sufficiently small to take [5\ = i in Eq. ([4]) (this 
only means that wrongly binned events are not shifted 
by more than one energy bin) and to expand the whole 
\ 2 function up to first order in b 2 , af, and bai. 

Combining the normalization and spectral uncertain- 
ties, Ti is given by 

Ti = (l + a l + b)[{N l+1 -N t )-b-(i + t Q + l) + N l ]. (5) 

Neglecting terms such as 0(b 3 ) and higher in \ 2 means 
neglecting terms of order 0(b 2 ) or G(ba,i) in Ti. With 
this approximation, our x 2 function becomes 



X 2 = E ( ]7 [(! + a * + + b ^ N ^ - N i> 



(i + t + ±)-Ni 



(6) 



This expression has to be minimized with respect to b: 
' >X = i \( a i + b ) N - + - W + *0 + \) 



db ^ Nj 

i 

■ [N i + (N i+1 -N i )(i + t + ^) 
b=- 



Si Ni T» 



where we have introduced the short hand notation 



7 i = N, + {N l+1 - Ni)(i + t + i). 



(7) 
(8) 

(9) 



There are two extreme cases: First, let us assume a very 
smooth energy spectrum, i.e. iVj+i — Ni <C Ni. Then, 
fa Ni, and it follows from Eq. (JSJ) that b is of the same 
order as a^. However, if the energy spectrum contains 
several large steps (as the Geo-neutrino spectrum does) , 
at least some of the differences N i+1 — Ni are sizeable, 
so ji 3> Ni. As 7, enters the denominator of Eq. ([5]) 
quadratically, while the numerator exhibits only a linear 
dependence, the fit value of b will be very small in this 
case (close to zero), i.e. the energy calibration uncertainty 
is essentially eliminated. 



3. APPLICATION TO EARTH MATTER 
EFFECTS ON SUPERNOVA NEUTRINOS 

On the one hand, the self-calibration effect is a nice 
tool to make an experimentalist's work more efficient. 
But to demonstrate that it can also help to do preci- 
sion measurements, we consider the matter effects of the 
Earth on supernova neutrinos as an example using the 
same LLSD as in Sec. [5J 

The fact that Earth matter has an effect on the 
spectrum of supernova neutrinos has been presented 
in 0, HU. A supernova core is essentially a neutrino 
blackbody source, so that the spectrum of F e 's produced 
in the supernova is thermal. However, it is important 
to keep in mind that the neutrinos that are produced as 
V e reach us as matter eigenstates V\, which means that 
the oscillations that appear inside the Earth are V\— V<i 
oscillations 16]. This has been taken into account in our 
simulations by modifying the source code of the GLoBES 
software accordingly. 

The effect of the oscillations in the Earth on the neu- 
trino energy spectrum is then, that one can see wiggles on 
the otherwise smooth spectrum. This can be shown more 
clearly by going at first to inverse energy units and then 
taking the squared modulus of the Fourier transform of 
the resulting spectrum, which gives the power spectrum, 
that points out the different oscillation modes. In the 
case of Earth matter effects, this results in one or more 
peaks, depending on whether the neutrinos traverse the 
core of the Earth or only the mantle. To extract these 
peaks from the experimental data, it is important to pre- 
cisely know their positions. This obviously requires an 
excellent energy calibration [l5l |. 

As example for our calculations, we have taken the 
"accretion-phase model I" from pj} • This gives the value 
a = 4.4 for the flux parameter as well as <i>(77 e )/$(77 x ) = 
0.8 (for details, see jl5|), where x stands for all other neu- 
trino flavours except T> e . For modelling the Earth matter, 
we have taken two different scenarios: a constant density 
equal to the average mantle density of 4.6 g/cm 3 accord- 
ing to the PREM profile [18] as well as a 3 layer approxi- 
mation with two times the mantle density, but a different 
core density of 11.8 g/cm 3 . In both scenarios, the neutri- 
nos travel through the whole Earth, which corresponds 
to a baseline of 12742 km. Note that our normalization 
is made for approximately 2000 events in the detector, 
which is a good example value according to [3] ■ 

This has been simulated with the GLoBES software 
package for different cases (cf. Fig. Q|. The first case 
is propagation in vacuum (black-dotted line). This does 
of course not have any effect on the power spectrum of 
neutrinos reaching the Earth from the supernova, simply 
because it is completey equivalent to a 12742 km longer 
propagation in space, which clearly makes no difference 
for a supernova whose expected distance from the Earth 
is several kpc. For propagation through matter, one can 
see - as expected - one or even more peaks, depending on 
the number of layers of constant density in the consid- 
ered Earth model. The red-dashed curve with a wrong 
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energy calibration (10% error) is clearly separated from 
the one with perfect calibration (green-solid). Of course, 
10% is much worse than this error would be in reality, 
but here just the Self-Calibration Effect is to be demon- 
strated which is more illustrative for a larger error. The 
crosses are data points coming from a simulation that 
has also been started with an initial calibration error of 
10%, but in that case, we have fitted the energy calibra- 
tion b to the Geo-neutrino background. This can be done 
since the SN neutrinos can be easily separated from all 
other events simply because of the narrow time window 
for neutrinos coming from a supernova. This \ 2 analysis 
then pulls the value of the energy calibration b down to 
zero because of the self-calibration effect. Then, shift- 
ing the energy of the events by the best-fit value of the 
calibration gives no difference to the case of perfect en- 
ergy calibration, which can be seen for both scenarios in 
Fig.Hl 

Note that in this calculation, we have assumed a per- 
fect extrapolation of the energy calibration to higher en- 
ergies. Geo-neutrinos have energies up to about 3.3 MeV 
(cf. Fig. [T]), while supernova neutrinos can have energies 
of several tens of MeV. Of course, in reality, this extrap- 
olation cannot work perfectly, but may be possible to a 
good enough precision. 

This shows that the Self-Calbration Effect can be a 
nice tool for measuring Earth matter effects on neutri- 
nos produced in supernovae. However, this is just one 
example, but of course, the effect applies to all situa- 
tions, where an accurate energy calibration is required 
and a background source with characteristic spectrum 
and good enough statistics is available. Possible applica- 
tions of the solar neutrino spectrum should also be con- 



sidered. 

4. CONCLUSIONS 

We have presented the interesting possibility for a 
Large Liquid Scintillator Detector in future neutrino ex- 
periments to use natural backgrounds with a character- 
istic spectral shape for the energy calibration. Due to 
their well-known spectrum, this could be Geo-neutrinos 
as source of i> e 's or solar neutrinos as source of i^'s. We 
have demonstrated that this effect is only tiny for re- 
actor neutrino experiments due to their high statistics 
and known spectral shape. However, we have also shown 
that the effect can be derived analytically and that it 
could be a nice tool in suitable scenarios, e.g. supernovae, 
where the correct energy calibration is crucial for de- 
tecting Earth matter effects on the neutrinos originating 
from such an explosion. The effect applies to all situa- 
tions with a suitable background source and high enough 
statistics. 
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Figure 1: The energy spectrum of Geo- neutrinos [13. fl3l] . Note that only the decay chains of uranium-238 and thorium-232 
produce neutrinos that are above the threshold of inverse f5 decay at 1.8 MeV. 



Influence of Geo-v's on # 12 (true value: sin 2 2# 12 =0.83) 
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Figure 2: Sensitivity of LENA to #12 for the assumed true value as a function of the baseline where the considered reactor has 
a thermal power of 0.5 GWth- The marked part is the region where the background self-calibration comes in. 
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Background Self-Calibration (SMALL scenario) 
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Figure 3: The rectangle from Fig. [2] drawn to a larger scale. Here, one can clearly see the seemingly paradoxial situation that 
a measurement with backgrounds can yield better results than one without. 
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Matter effects on Supernova neutrinos 



Constant density approximation 
AL=2R Earth 



Matter effects on Supernova neutrinos 
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Figure 4: Matter effects on the power spectrum of supernova neutrinos for a constant density profile and a 3-layer approximation. 
The neutrinos always propagate through the whole diameter of the Earth. Considered are the cases of propagation through 
vacuum and through matter, while the latter is subdivided in the case with perfect energy calibration and the one with a 
10% calibration error. Furthermore, for the wrong calibration, we have in one case performed a x analysis and shifted the 
events with the fittet value of the energy calibration. Due to the self-calibration effect, this analysis has again given the correct 
adjustment of the energy, which causes the corrected data points to lie exactly on the line of perfect energy calibration. 



